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ABSTRACT 

We report on the first XMM-Newton observation of tlie Vela-like pulsar 
PSR B2334+61. Spectral analysis reveals soft X-ray emission, with the bulk 
of the photons emitted at energies below ~ 1.5 keV. We find that the spectrum 
has a thermal origin and is well-fitted with either a blackbody or a magnetized, 
pure H atmospheric model. In the latter case, for a neutron star with a radius of 
13 km and a magnetic field of 10^^ G, the best-fit gives an hydrogen column den- 
sity Nh = 0.33 X 10^^ cm~^ and an effective temperature T^^ = 0.65 x 10^ K, as 
measured at Earth. A comparison of the surface temperature of PSR B2334+61 
obtained from this fit with cooling curves favor a medium mass neutron star with 
M ~ I.45M0 or M ~ I.6M0, depending on two different models of proton su- 
perfluidity in the interior. We do not detect any pulsed emission from the source, 
and determine an upper limit of 5% for the modulation amplitude of the emission 
on the pulsar's radio frequency. 

Subject headings: pulsars: individual (PSR B2334-f61) — stars: neutron — X- 
rays: stars 



1. Introduction 

Young pulsars, neutron stars in supernova remnants (SNRs) and cooling neutron stars 
make up a subset of rotation-powered radio pulsars that can be observable at shorter wave- 
lengths. Measurements of the high energy radiation are crucial as they provide important 
insights into the emission processes and the neutron star physics. 
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Through observation and theoretical modehng it is generally accepted that a combi- 
nation of emission mechanisms are responsible for the detected X-ray flux from rotation- 
powered pulsars (for a review see e.g. Becker & Pavlov 2001; Becker & Aschenbach 2002). 
Soft thermal radiation is produced by cooling of the surface of the neutron star (cf. Green- 
stein & Hartke 1983; Romani 1987; Pavlov et al. 1995), while the acceleration of particles 
in the neutron star magnetosphere generates non-thermal radiation (cf. Michel 1991; Be- 
skin, Gurevich & Istomin 1993, and references therein). A harder thermal component can 
be produced by hot spots at the polar regions (Kundt & Schaaf 1993; Zavlin, Shibanov & 
Pavlov 1995), and the presence of a synchrotron nebula may contribute another source of 
non-thermal radiation (cf. Arons & Tavani 1993). 

The dominant emission mechanism is related to the age of the pulsar (see e.g. Becker et 
al. 1996; Pavlov et al. 2002). In pulsars older than 10^ yr standard coohng scenarios predict 
that the temperature of the neutron star surface will be too low to generate detectable 
thermal X-ray emission (Nomoto & Tsuruta 1987). However, the source may be observable 
if reheating via rotational energy loss or heating of the polar caps takes place (Brinkman, 
Ogelman & Aschenbach 1985). The middle-aged pulsars (10"^ — 10^ yr) exhibit spectra that 
can be described by thermal emission from the surface of the neutron star. In some cases a 
power-law component is also detected (e.g. Becker & Aschenbach 2002; Pavlov et al. 2002). 
In pulsars younger than 10^ yr the strong magnetospheric emission prevails over the weaker 
thermal radiation. By comparing measurements of the soft X-ray flux emitted from the 
neutron star surface with thermal evolution calculations it is possible to investigate the 
physical processes that occur in the interiors of these objects (Schaab et al. 1999). 

Although many rotation-powered pulsars have now been observed in X-rays, when one 
considers the number of known radio pulsars, the sample of sources detected at X-ray energies 
is still small. The prospects of detecting a greater number of these objects is however 
improving with missions like XMM-Newton and Chandra. It is now possible, even for faint 
and distant objects like the one presented in this paper, to discriminate whether the dominant 
mechanism is thermal or non-thermal which in turn permits the first estimates of temperature 
and/or power law index. In particular, the detection of emission due to cooling is essential 
to determine where the pulsar lies on the thermal/evolutionary diagram, with the main 
objective being able to constrain observationally the equation of state of matter at supra- 
nuclear densities. Gurrently, there are very few sources for which this is possible, therefore 
any new source that can be added to the parameter space is important. 

One such object that should fall into the category of pulsars for which the cooling 
neutron star is the dominant source of emission is PSR B2334-I-61. This Vela-like pulsar is 
located at a distance oi D — 3.ll5;o kpc (Gordes & Lazio 2002). Analysis by Kulkarni et al. 
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(1993) implies that the pulsar is associated with the SNR G114.3, making the source one of 
the oldest (spin-down age of ~ 4 x 10^ yr) that is still linked to a SNR. Hence, PSR B2334+61 
is a key object to study with relation to the thermal/evolutionary parameter space. The 
spin-down rate (P = 191 x 10~^^ s s"^) of the 495 ms pulsar indicates that the magnetic field 
of PSR B2334+61 is ~ lO^^ q pgj^ B2334+61 was originally detected in a short (8 ks) 
i? 05*^ r pointing (Becker et al. 1993, 1996) but the low statistics prevented a comprehensive 
investigation of the spectral properties of the source or a meaningful temporal analysis. Here 
we report on the first XMM-Newton observation of PSR B2334-F61. 



2. Observations and data reduction 

PSR B2334+61 was observed with XMM-Newton on 2004 February 12. For the spectral 
and timing analysis we used data from the European Photon Imaging Camera (EPIC) in- 
struments: the EPIC MOS detector (Turner et al. 2001) and the EPIC-PN detector (Striider 
et al. 2001). Both MOS instruments and the PN were configured in full frame mode and we 
used the thin filter. The MOSl and M0S2 observations had total exposure times of 47.5 
ks. The EPIC-PN buffer was filled during the observation resulting in an exposure time of 
45.8 ks. Wc reduced the EPIC data with the XMM-Newton Science Analysis System (SAS 
version 6.0.0). 

In order to maximize the signal-to-noise ratio for our XMM-Newton observation, we 
filtered the data to include only single, double, triple and quadruple photon events for the 
MOS, and only single and double photon events for the PN. We included photons with 
energies only in the range 0.3 — 10 keV. 

To verify that there is no contribution to the measured emission from the remnant, 
wc compared the radial profile of the pulsar emission with the XMM-Newton point-spread 
function for EPIC-PN at 1.5 keV, which we generated using the King profile parameters 
included in the XMM-Newton calibration file "XRT3_XPSF_0006.CCF.ph" ^ Wc find that 
the emission we detect from PSR B2334-I-61 is consistent with that from a point source. 

Data have been extracted using a circular region of radius 30", centered on the radio 
pulsar's position. Wc obtained 600 counts from MOSl, 616 from M0S2, and 2384 from 
PN. By using a region of equal size offset from the pulsar's position, we found that the 
background contributes 439 counts in MOSl, 449 in M0S2, and 1846 in PN. We detect 
X-rays from PSR B2334+61 at a level of 5a in MOSl, 5.1a in M0S2, and 8.3a in PN. The 
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extraction radius we used encircles 85% of the energy from the source in MOS and 83% in 
PN. The measured fluxes reported in Table 1 have been corrected accordingly. 

3. Spectral Analysis 

In order to perform the spectral analysis, we extracted the pulsar spectrum from the 
MOSl, M0S2 and PN event files. Data were filtered to exclude events that may be incorrect, 
for example those next to the edges of the CCDs and next to bad pixels. The three spectra 
were regrouped by requiring at least 50 counts per spectral bin. For each spectrum we created 
a photon redistribution matrix (RMF) and ancillary region file (ARE). The subsequent 
spectral fitting and analysis was performed using XSPEC, version 11.3.1. 

Since essentially no emission is detected at energies > 1.5 keV, we modeled the MOS 
spectra over 0.5 ^1.5 kcV and the PN spectrum over 0.3 — 1.5 keV (taking into account the 
calibration uncertainties for the MOS at low energies). We fitted the combined MOS and 
PN spectra using different models: initially a power-law, a blackbody and a combination of 
both, each modified by photoelectric absorption (wabs in XSPEC). The best-fit results are 
summarized in Table 1. The values indicate that the spectrum is better represented by 
a thermal (blackbody) model, rather than by an absorbed power-law. The additional power 
law component in the third (and fifth) fit is poorly constrained, and not required: by applying 
an F-test, we find an F-statistic of 1.6 and probability of 0.2. The best-fit blackbody fit gives 
Nh = 0.26 X 10^2 cm-2 and T°° = 1.62 x 10^ K (where T~ is the blackbody temperature 
measured ad Earth^; see Figure 1, first and second panels). We notice that, consistently, this 
value for the hydrogen column density is lower than the Galactic one in the source direction, 
which is Nh ~ 0.84 x 10^^ cm"^, while the same argument can only be applied marginally 
(within the lower limit of the error bars) to the value of Nh inferred from the power-law 
fit. The radius of the emitting region, as derived from the blackbody fit using the Cordes & 
Lazio (2002) distance, is Rem = l-66l|]:i km. 

In order to attempt a different representation of the thermal emission detected from 
PSR B2334-I-61, we used a magnetized, pure H atmospheric model (nsa in XSPEC, for details 
see Pavlov et al. 1995), by fixing the magnetic field aX B — 10^^ G, the value inferred from the 
radio timing measurements of the source. The neutron star mass and radius (measured at 
the source) were fixed at M^s — 1.4M© and it! = 10 km, respectively, although the spectral 



^We note that can be obtained by red-shifting at Earth the effective temperature measured at 
the neutron star surface, Tg//, through the relation = Te//(1 — 2.952M„s/i?)°'^. T^ff is one of the 
best-fitting parameters of the NSA model. 
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fit is not particularly sensitive to these parameters (see §5 for an a posteriori estimate of 
the mass). We found that this model also provides an excellent representation of the pulsar 
spectrum, and the best-fit parameters are Nh = 0.42 x 10^^ cm~^ and = 0.58 x 10® K, 
where T^j- is the effective temperature as measured at Earth^ (see Table 1). However, this 
fit gives a pulsar distance of = 1.1 ±0.6 kpc, lower than that inferred from the value of the 
electron density {D = 3.11?;q kpc, Cordes & Lazio 2002). Although the pulsar distance and 
radius are poorly constrained by our data, we find that a better agreement can be reached 
by using the same atmospheric model and fixing the neutron star radius at a slightly larger 
value of 13 km. In this case, we find Nh = 0.33 x 10^^ cm-^ and T^f = 0.65 x 10® K 
(see Table 1 and Figure 1, third and fourth panels). The distance resulting from this fit 
is D = 3.2 ± 1.7 kpc, which is consistent with the electron density distance. Finally, we 
verified that even when the thermal component is modeled by a neutron star atmosphere, 
an additional power law is not required (last fit in Table 1). 

We note that the residuals from the model fits shown in Figure 1 are relatively large 
at ~ 0.6 — 0.7 keV. To try and reduce the residuals we fitted the data using the previous 
models, each modified by an absorption, namely phabs or tbabs in XSPEC. We find that 
the resulting fits are not improved with respect to our previous results. 

4. Timing Analysis 

For the timing analysis we only used FN data, extracted by applying the filtering criteria 
and extraction region previously described in §2; the filtered file was then barycentrically 
corrected. In order to search for an X-ray modulation at the PSR B2334+61 spin period, 
we first determined a predicted pulse period at the epoch of our XMM-Newton observations, 
assuming a linear spin-down rate and using the radio measurements (Dewey et al. 1985; 
Hobbs et al. 2004). We calculate P = 0.495355469 s (/ = 2.0187523 Hz) at the midpoint of 
our observation (MJD 53,047.6). As glitches and/or deviations from a linear spin-down may 
alter the period evolution, we then searched for a pulsed signal over a wider frequency range 
centered on f — 2.01875 Hz. We searched for pulsed emission using two methods. In the 
first method we implement the test (Buccheri et al. 1983), with the number of harmonics 
n being varied from 1 to 5. In the second method we calculate the Rayleigh statistic (de 
Jager 1991; Mardia 1972) and then calculate the maximum likelihood periodogram (MLP, 
Zanc et al. 2002) using the C statistic (Cash 1979) to determine significant periodicities in 
the data sets. 

We do not find any significant peak near to the predicted frequency with either method, 
either using the whole 0.3 — 10 keV energy band or restricting the search to the 0.3 — 1.5 keV 
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band. By folding the light curve of PSR B2334+61 on the radio frequency and fitting it with 
a sinusoid, wc determine an upper limit for the pulsation of 5% in modulation amplitude 
(defined as {Fmax-Fmin)l {Fmax+Fmin) where F^ax and F^in are the maximum and minimum 
of the pulse light curve) . 

5. Discussion 

We have presented the results from the first XMM-Newton observation of PSR B2334+61. 
The source has been positively detected in all EPIC instruments, although the X-ray emis- 
sion is very faint and the spectrum does not have a high enough signal-to- noise to make 
a detailed multi-component fit. However, single component fits allowed us to discriminate 
between a thermal (blackbody-like) or non-thermal (power law-like) nature of the domi- 
nant emission mechanism. We find that the spectrum is well represented by either thermal 
blackbody-like emission from a small emitting area (e.g. a hot polar cap) of size ~ 2 km or 
emission from a pure H atmosphere of a neutron star with i? ~ 13 km and magnetic field 
B ~ 10^3 Q The EPIC flux, measured in the 0.3 - 10 keV band, is (9.21°;^) x 10"^^ erg 
cm -2 s-i and (l.rtil) X 10"^^ 

erg cm ^ s ^, in the two cases respectively. Both values are 
slightly higher than that previously inferred by Becker et al. (1996), (7.1 ± 0.2) x 10~^^ erg 
cm~^ s~^, although we note that due to their limited counts ROSAT data could only be 
analyzed by making an assumption a priori on the spectral shape. In both cases, the value 
of Nh is consistently lower than the total galactic one in the pulsar direction. A fit with a 
power-law model alone is statistically worse, and we find little evidence for the presence of a 
non-thermal component in the spectrum in addition to the thermal one. The upper limit for 
the flux contribution from the power- law in the energy range 0.3 — 1.5 keV to the blackbody 
flt is 3%, and to the magnetized atmospheric model flt is << 1%. We do not detect X-ray 
pulsations corresponding to the radio signal, to a limit of 5% in modulation amplitude. 

Although both thermal flts are conceivable options, we tend to prefer the atmospheric 
model representation based on physical grounds. In this case, magnetic effects for a fleld 
strength of order ~ 10^^ G (which is the value inferred from measurements of the source spin 
period and period derivative in the radio band) are consistently accounted for in the radiative 
transfer computation. Moreover, a value for the radius of ~ 10 — 13 km is in agreement with 
the prediction of several neutron star equations of state (Lattimer & Prakash 2001) and we 
found that this parameter can be successfully adjusted to make the distance inferred from 
the spectral fit consistent with that obtained from the electron density value of this pulsar. 

To date, thermal emission has been detected in only very few radio pulsars: PSR B0656+14 
(Possenti et al. 1996), PSR B1055-52 (Pavlov et al. 2002), PSR J0437-4715 (Zavhn et al. 
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2002), PSR J0538+2817 (McGowan et al. 2003), Geminga (Halpern & Wang 1997), Vela 
(Pavlov et al. 2001), and PSR B1706-44 (Gotthelf et al. 2002, McGowan ct al. 2004). In 
the case of the first three objects the thermal emission detected above ~ 0.5 keV is more 
likely to originate from a hot-polar cap. These sources are in fact old pulsars, at a more 
advanced stage of their cooling history and their surface emission should peak at much lower 
(UV) energies. This idea was strengthened by the detection of a further thermal component 
below 0.7 keV (Pavlov et al. 2002) in the spectrum of the brightest of them, PSR B0656+14. 
The Vela pulsar and PSR B 1706-44 are younger (r ~ 10^ yrs) and are the only radio active 
sources for which the thermal component observed in the soft X-rays is well explained by a 
magnetized cooling atmosphere (Pavlov et al. 2001, McGowan et al. 2004). When this model 
is assumed instead of a blackbody, the inferred radius is in agreement with a neutron star 
equation of state. The only other neutron stars whose thermal component is better described 
by an atmospheric model, and for which this interpretation resolves all the inconsistencies 
which follow from the blackbody interpretation, are the radio-silent neutron stars IE 1207-52 
(Zavlin et al. 1998) and RX J0822-4300 (Zavlin et al. 1999). 

The XMM-Newton observation reported here allows us to add another entry to the 
list. If indeed the X-ray emission detected from PSR B2334+61 is originating in the cooling 
atmosphere of the neutron star, our estimate of the effective temperature allows us to localize 
the object in the neutron stars thermal evolutionary diagram. 

Our knowledge of neutron star interiors is still uncertain and accurate measurements 
of the neutron star surface temperature arc particularly important to constrain the cooling 
models and provide information on the physics of the neutron star. Roughly speaking, theo- 
retical models predict a two-fold behavior of the cooling curves depending on the star mass. 
In low-mass neutron stars neutrino emission is mainly due to a modified Urea process and 
nucleon-nucleon bremsstrahlung. These are relatively weak mechanisms and produce slow 
cooling. In stars of higher mass the neutrino emission is enhanced by a direct Urea process 
(or other mechanisms in exotic matter), therefore these stars cool down much faster {fast 
cooling regime). To date (see Yakovlev et al. 2002 for a discussion) it has been realized that 
simple models which do not account for proton and neutron superfluidity fail in explain- 
ing the surface temperatures observed in many sources, unless objects such as e.g. Vela, 
Geminga, RX J1856-3754 do have exactly the critical mass that bounds the transition be- 
tween the very different slow cooling and fast cooling regimes. This unlikely fine-tuning is not 
required if the effects of nucleon superfluidity are accounted for. In particular, models with 
proton superfluidity included predict an intermediate region between fast cooling and slow 
cooling curves, which is expected to be populated by medium mass neutron stars (roughly 
with M between 1.4 and 1.65 Mq). Although the full picture only holds if, at the same 
time, neutron superfluidity is assumed to be rather weak, it is still interesting that many 
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neutron stars (as IE 1207-52, Vela, RX J1856-3754, PSR 0656+14) have a surface temper- 
ature which falls in such a transition region (Yakovlev et al. 2002). In turn, this means 
that measuring the surface temperature allows us to "weigh" neutron stars (Kaminker et 
al. 2001). As we can see from the first two panels of Figure 2 in Yakovlev et al. (2002), 
assuming an age of logr = 4.6, the surface temperature of PSR B2334+61 derived from 
the blackbody fit is even higher than the upper cooling curves i.e. those corresponding to 
the slow cooling regime. However, the surface temperature \ogT°° — 5.9 obtained by fitting 
with the magnetized model and it! = 13 km falls well within the above mentioned transition 
region of medium mass neutron stars. The mass of PSR B2334+61 should then be ~ 1.45M0 
or ~ I.6M0, depending on the kind of proton superfluidity assumed in the model (Ip and 
2p respectively). Indeed, the measured properties of PSR B2334+61 reported here are in 
remarkable agreement with that measured for PSR B1706-44 (McGowan et al. 2004). 

As mentioned in the introduction, we are now in the position to make the first experi- 
mental studies of surface temperatures of isolated neutron stars. However, although most of 
our insight of neutron star temperatures and interior rely on them, at present these studies 
have to be considered as pioneering. Only a few sources are currently available for this kind 
of exercise, so every newly discovered candidate is important. At present. X-ray spectra 
of thermal emission from NSs are fitted with either a blackbody or magnetic atmosphere 
models. Although the latter definitely represent a substantial improvement, inasmuch they 
include most of the relevant physics, a large amount of work remains to be done before they 
can be claimed to be fully satisfactory. The same is true for the coohng curves: progress 
has been made in improving these models but they can not yet be considered as completely 
realistic. Besides other effects, we make the caveat that both spectra and cooling curves are 
computed assuming 1-dimensional transfer of radiation/heat in a single effective tempera- 
ture, a single magnetic field zone and neglect all effects of more realistic neutron star thermal 
and magnetic surface distributions (sec Zanc & Turolla 2005, Page ct al. 2004, Blaschkc ct 
al. 2004). However, what values to include for these parameters is not necessarily obvious as 
it is not completely clear what the NS surface temperature distribution is even in the case of 
a simple dipolar magnetic field (see Geppert, Kiiker & Page 2004). We also note that inter- 
preting the temperatures obtained from the spectral fits in the context of theoretical cooling 
curves relies on the true age of the pulsar being the same as the characteristic spin-down 
age, which may not be valid. Our results for PSR B2334-I-61 provide information that helps 
to constrain the current models and will enable more realistic models to be produced in the 
future. 

This work is based on observations obtained with XMM-Newton, an ESA science mission 
with instruments and contributions directly funded by ESA Member States and NASA. 
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Fig. 1. — XMM-Newton PN and MOS spectra of the X-ray emission detected from 
PSR B2334+61. First panel: data (crosses) and best-fit blackbody model (solid line) for 
the parameters given in Table 1. Second panel: Difference between the data and the black- 
body model. Third panel: data (crosses) and best-fit magnetized, pure H, atmospheric 
model (solid line) for the parameters given in Table 1 (NSA^). Fourth panel: Difference 
between the data and the atmospheric model. In the electronic edition the PN, MOSl and 
M0S2 data points, and best fit model where appropriate, are colored blue, red and green, 
respectively. 
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Table 1. Spectral fits to the X-ray emission from PSR B2334-I-61 
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PL 


0.901°:°^ 






38/25 = 


1.5 


{<oAtf^^) X 10- 


-11 


BB 


0.261°;^^ 




-'-■"^-0.23 


27/25 = 


1.1 


(9.21°:^) X 10- 


-14 


BB+PL 


0.431°;?^ 


r, r,+3.0 
^■^-1.4 


-1 97+0.35 
-0.58 


24/23 = 


1.0 


(3.1^}i) X 10- 


-13 


NSA'^ 


^•"^^-0.04 




U.05_o.25 


28/25 = 


1.1 


(3.4^1°) X 10- 


-13 




'-'•'J'-'-O.lO 






31/25 = 


1.2 


(1.7i;i) X 10- 


-13 


NSAi^+PL 


•-"•^".Q Q4 


q 4+0.6 


n 76+0-09 

'-'•'"-0.65 


34/23 = 


1.5 


(1.0l?:S) X 10- 


-13 



Note. — Uncertainties are all at 90% confidence. Models used: "PL" indicates 
an absorbed power law with photon index F; "BB" indicates an absorbed black- 
body emitting at a temperature T; "NSA" indicates an absorbed, magnetized, 
pure H atmospheric model with B — 10^^ G and effective temperature as 
measured at Earth (Pavlov et al. 1995). For the NSA model we fixed the neutron 
star mass at Mns = 1.4M0, and the radius at 10 km (NSA*") or 13 km (NSA^). 
The last column is the unabsorbed flux in the 0.3 — 10 keV band. 



